Exact Ground State and Finite Size Scaling in a 
Supersymmetric Lattice Model 
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We study a model of strongly correlated fermions in one dimension with extended N=2 super- 
symmetry. The model is related to the spin S — 1/2 XXZ Heisenberg chain at anisotropy A = — 1/2 
with a real magnetic field on the boundary. We exploit the combinatorial properties of the ground 
state to determine its exact wave function on finite lattices with up to 30 sites. We compute several 
correlation functions of the fermionic and spin fields. We discuss the continuum limit by construct- 
ing lattice observables with well defined finite size scaling behavior. For the fermionic model with 
periodic boundary conditions we give the emptiness formation probability in closed form. 
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Supersymmctry is well motivated in high energy 
physics where it offers a partial solution to fine-tuning 
problems and improve gauge coupling unification [j]. It 
also appears in condensed matter models although at 
a less fundamental level. Examples are disordered sys- 
tems Q and models of strongly correlated electrons like 
extended Hubbard ,3j or t-J models Q where supersym- 
metry relates fermionic and bosonic composite operators. 

A typical consequence of unbroken supersymmetry is 
the prediction of the ground state energy. This is not 
sufficient to compute the ground state wave function, the 
relevant quantity for the calculation of vacuum expecta- 
tion values. It is quite natural to ask whether supersym- 
metry and the knowledge of the ground state energy are 
useful to this purpose. 

In this Letter, we analyze the problem in a recently 
proposed one dimensional model of itinerant fermions p| 
with two supercharges obeying with the Hamiltonian an 
extended N=2 supersymmetry algebra. The knowledge 
of the ground state on large finite lattices is important 
to study its continuum limit where the model is expected 
to describe a minimal superconformal series. 

We shall be interested in boundary effects and thus 
consider mainly free boundary conditions, see p for 
the periodic case. The model can be mapped to the 
integrable open XXZ Heisenberg spin 1/2 chain with 
anisotropy A = —1/2 and a suitable real surface mag- 
netic field. In principle, Bethe Ansatz techniques could 
be applied. The supersymmetry inherited from the 
fermionic model should allow to compute Baxter function 
whose zeros give the Bethe quantum numbers How- 
ever, the procedure is rather involved with open bound- 
ary conditions as discussed in p}- 

Here, we pursue a different approach starting from the 
following remarks. The XXZ chain at A = —1/2 is inte- 
grable for a large class of boundary conditions. In some 
specific cases (e.g. twisted or U q (sl(2)) symmetric ones) 
several remarkable conjectures have been claimed about 
the combinatorial properties of the ground state wave 
function 0, E| . They arise from the relation between the 



XXZ chain and Temperley-Lieb loop models Ell- 
in this Letter, we first show that similar features are 
present in the fermionic model and in the related XXZ 
chain with surface magnetic field. We then explain how 
non standard number theoretical methods can be used 
to obtain exact expressions for the ground state wave 
function of the fermionic model on long chains. Finally, 
we analyze in some details the physical properties and 
finite size scaling (FSS) behavior of the ground state. In 
particular, we discuss the continuum limit and propose 
a way to extract scaling fields from the fermionic model 
and the associated XXZ Heisenberg chain. 

The model is defined on a one dimensional lattice 
with L sites and free boundary conditions. Let Cj, c\ be 
spinless fermionic creation annihilation operators with 
algebra {c i; c}} = <Sy, {ci,Cj} = {c],ct} = 0. We de- 
note by Ni the set of nearest neighbours of site i. The 
projector over states with a hard core condition forbid- 
ding occupancy around site i is "Pi = YijeN (1 ~ n j)> 



with n, 



c]c 3 . 



vj — ^jisj. Let Q + be the supersymmetry charge 
Q+ = Ei4 Vi. and Q~ = (Q+)t. The operators Q ± 
are nilpotent. The Hamiltonian 



H = {Q+,Q-}=J2J2 V > c lcj'P J 



(i) 



is by construction symmetric. We restrict ourselves 
on the subspace Hl,f of states with F fermions and no 
adjacent occupied sites which is an invariant subspace of 
the full Fock space under the action of . We work 
in the basis of simultaneous eigenstates of the number 
operators rij. The structure of eigenstates of H follow 
from supersymmetry. The energy is non negative and all 
energy eigenstates with E > are doublets connected 
by the action of Q ± . The zero energy states are singlets 
annihilated by Q^. They are supersymmetric ground 
states. A cohomological analysis 5] shows that there is a 
unique zero energy state for L mod 3 = 0,2 and none for 
L mod 3 = 1. In the following we shall consider the case 
L = 3n. The unique ground state has then fermion num- 
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ber F = n. The dimension of Wj,,£/3 increases rapidly as 
L — > oo. It should be clear that the determination of the 
unique ground state |^ ) is non trivial on large lattices. 

At low L, the ground state can easily be obtained. In- 
spection of the explicit wave functions reveals a remark- 
able fact. Indeed, the ground state can always be written 
in the very special form 



1010- ■ -1000- ■ -0) +^x s |s) 



(2) 



2F terms F terms 



where x s € Z and s runs over all states (diagonal in 
the occupation number basis) with the exception of the 
single state explicitly written. Lattice parity halves the 
number of independent coefficients. We used it as a con- 
sistency check. The integrality property of {x s } is def- 
initely a non trivial Ansatz on the ground state wave 
function, and it would be elusive in a Bethe Ansatz ap- 
proach. Similar results already appeared in the literature 
for the XXZ Heisenberg model atA = -l/2@,@ that 
is indeed closely related to the present fermionic model. 

Let us assume the integrality property (0) as a working 
hypothesis. The equation H\^ ) = reduces to a linear 
problem of the special kind 



Ax = b, AeZ dxd ,beZ 



vlxd 



(3) 



where x is the vector of unknown coefficients x s and 
d = dim7YL.F — 1- We know that this problem admits 
an integer solution x £ Z lxd and it seems reasonable to 
be able to find it exactly with modest effort. This is pos- 
sible due to a well known technique in cryptology. The 
problem can be solved by working in the finite field Z p of 
integers modulo a prime p. To solve Eq. ©, we choose 
a large prime p and first determine a solution modulo p 
by the Lanczos algorithm over finite fields ■ We build 
the sequences {&i}o<i<<i and {ci}i<i<d where the initial 
values are 



b Q = b, ci = Ab, h = Cl - — ^b, (4) 

o • Ci 

and the sequences are generated by iterating 



c l+1 = Abi, 

2 

t c i+l , c i ' c i+l ■ 

b l+ i = c i+ i - bi / 

bi ■ c l+ i c t ■ Oj_i 

The solution is finally obtained as 
£5 h-b , 



bi ■ Ci + i 



(5) 
(6) 

(7) 



In Eqs. I|4I7I) . all arithmetic operations, in particular di- 
visions, are done in Z p . The above algorithm exploits 
the sparsity of A, and is not memory expensive since 
requires a storage 0(d). After a certain number of solu- 
tions modulo (large) primes {pi} have been found, they 



are combined together to give the exact solution. This is 
done by applying the Chinese Remainder Theorem [l^ . 
By the above procedure, we have determined the exact 
ground state for L = 3n < 30 (dim ^30,10 = 352716). 
The complete expressions (i.e. the sequences {x s }) are 
available on request. 

We now illustrate the main physical properties of 
the results. For any observable O we denote (O) = 
(^olCl^o)- Several remarkable combinatorial features 
appear immediately leading to conjectures in the spirit 
of 0. We find that the dominant state in the ground 
state is always |r) = |010010010 ■ ■ • ) with coefficient 
max|a; s | = X\ T ) = A 8 (2L/3 + 2) where iV 8 (2n) = 

llfc^iC 3 *; + l)(2/c)!(6fc)!/((4/c)!(4fc + 1)!) is the number 
of cyclically symmetric transpose complement plane par- 
titions 0. Also, the squared norm in our normalization 
is 1+J2 S x 1 = A r s(2L/3+2)Ay(2 J L/3+3), where A v (2n+ 
1) = 2" n n^ =1 (6A;-2)!(2A;-l)!/((4fc-l)!(4fc-2)!)is the 
number of vertically symmetric alternating sign matrices 
We checked these expressions for all L < 30. Similar con- 
jectures can be claimed on correlation functions, as in the 
recent work 0] . A simple example is the expectation of 
the potential energy Vi). The operator (1/L) ^\ T^i 
is diagonal on states |n) and is minimum on |r) where 
it attains the value F/L — 1/3. Our data are consistent 
with the simple formula 

1 5L + 21 

3 4TTT5- (8) 

1 

The asymptotic value is 5/12, slightly larger than 1/3, 
due to the subdominant states |n) ^ |t) appearing in 
I 'So)- Another example is the fermion density on the 
boundary (rti). We find 

L(10L + 33) 



"1 



2(4L + 9)(4 J L+ 15)' 



(9) 



The asymptotic value 5/16 is smaller than 1/3, a plausi- 
ble fact since the dominant state |r) has no fermions on 
the lattice boundary. We now examine the FSS behavior 
of simple correlation functions. In Fig. we show the 
expectation value (rife) — 1/3 on the L=30 lattice. There 
is a clear Z 3 substructure similar to that observed in the 
phase diagram of bosonic models with hard-core repul- 
sion |l4|. The average number of fermions in the three 
branches can be defined as Fk — J2a mod 3=fe( n i)- We 
have F = Fi and F + Fi + F 2 = F. Asymptotically, 
at large F = £/3, the upper branch has the asymptotic 
value F 2 = \{F + 1) + 0(F~ l ). The three sublattices 
containing sites k with fixed k mod 3 appear to recon- 
struct a smooth curve for (n^) as L increases. The Z3 
structure can be exploited to construct scaling fields in 
the continuum limit. Let us define the effective length 
L = L/3 + 1, and set k± = (L ± l)/2. Our data for (nu) 
suggest to test the following FSS laws 

(n k ) - 1/3 = / + f (k - k + )/t] I'", k mod 3 = 2, (10) 
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lation function Gij — (niTij) — (n,-)(nj). On the left side 



k mod 3 = 1 
a-a k mod 3 = 2 

k mod 3 = 



FIG. 1: Z3 structure of the expectation (rife). Dashed lines 
connect branches with the same value of k mod 3. 
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FIG. 3: Boundary and bulk pair correlation function at L=30. 
Algebraic decay of the k mod 3 = 1 component of the bound- 
ary correlation. 
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FIG. 2: FSS of the occupation number (n/c). 



(n fe +n fe+2 )-2/3=/_ y(k — &_)/LJ L~ v , k mod 3=1, 

where v and v' are unknown exponents. The best 
collapse of data at different L is obtained with v = 
v' = 0.33(2) and is shown in Fig. J2J. Notice that for 
k mod 3 = 1, the separate (n^) and (rifc+2) form a parity 
doublet. They should be related in the continuum limit 
to the left and right moving parts of a single field. The 
value of v, v' is close to 1/3. In principle, this information 
is useful in the effort of identifying the proposed scaling 
fields with the operator content of candidate supercon- 
formal field theories describing the continuum limit of 
the model. The next simple observable built with local 
fermionic fields is the density-density (connected) corre- 



of Fig. Q we show the boundary correlation function 
Gi^k on the L=30 lattice. Again the Z3 structure is evi- 
dent. On the right side, we fit G\ t k on lattices of various 
sizes with the simple form ag+ai/k with good agreement. 
The term ao is a small size effect decreasing with L. A 
detailed scaling analysis shall be reported elsewhere. 

The fermionic model can be mapped quasi-locally to 
the XXZ spin chain at anisotropy A = —1/2 with a spe- 
cific real surface magnetic field 0|. The structure of 
Bethe equations in the two models is closely related, but 
the fermionic and spin correlation functions can be quite 
different. Following |l5| . we map a fermionic configu- 
ration {n} to a well defined spin one {a} with the di- 
chotomic variable = ±1. The map is such that the 
spin chain has length 2L/3 + 1. The expectation value of 
the local spin is shown on the left side of Fig. QJ. The 
structure is now Z3 — > Z2 and (ofc) can be split into two 
independent channels with k mod 2 = 0,1. Eqs. (|8I9|) give 
the boundary spin (cti) = 2{n\) — 1 and the numbers of 
pairs of spins up or down 



(iV TT ) 



L(L - 3) 
3(4L + 9)' 



{Nil) 



L(L + 6) 
3(4L + 15)' 



(11) 



Again, it is possible to separate the two components in 
order to build well defined scaling fields. We test 



(a k )=g + Uk-L)/L\ 



L 



-0 



{o-k) = g- 



((k-L)/L) 



L 



k mod 2 = 0, 



k mod 2 = 1. 



(12) 



(13) 



The best values of the exponents are (3 = 0.66(1), f3' = 
0.80(1). The accuracy of the corresponding FSS laws 
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FIG. 4: (left) Z 2 structure, (right) FSS of (a k ). 
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-fe+1 



2/3 c f-v/3/2 



-5/36 



(17) 



where the constant c takes the form 



exp 



5e" 



36 



sinh(5t/12)sinh(t/12)\ dt 

- J — ( 18 ) 



sinh^/2) 



t 



To conclude, we have shown that the knowledge of the 
ground state |^o) on moderately large finite lattices al- 
lows a precise FSS analysis of correlation functions guid- 
ing the detailed identification of the continuum limit of 
the model (QJ. The exact |^o) is also an effective heuristic 
tool to derive closed expressions for correlation functions. 
We acknowledge conversations with P. Fendley and J. De 
Boer. 



is shown on the right side of Fig. Notice that (3, 

[3' are quite close to the simple rationals 2/3 and 4/5. 
This suggests again possible simple identification of these 
scaling fields with superconformal operators. 

We want to stress at this point that the same methods 
can be applied to any model with integrality properties 
like (J2J). In particular, a suitable version of J2J holds 
in the model Q with periodic boundary conditions. To 
give an example of the efficacy of our number theoretical 
techniques in this context, we have computed the unique 
supersymmetric ground state for the model with periodic 
boundary conditions and L — 3F + 1, up to L = 28. A 
relevant non-trivial observable is the emptyness forma- 
tion probability E k = (IIi=i(l ~ n i))- The first values 
are fixed from translation invariance and the constraint 
niUi+i — which defines H.l,f 



Ex 



2F+1 



F + l 



3F + 1' 3F+1 
For k > 2 we checked the validity of the relation 



Ek-i 



(A -2)1(3* -5)! 



F + v 



(2fc-3)!(2fc-4)! v ±*-_ k 2F + v' 



fc-i 

n 



(14) 



(15) 



We guessed the very specific form of Eq. (|15[) from similar 
conjectured relations proposed for the finite XXZ model 
with twisted boundary conditions |8|]. The thermody- 
namical limit F — > oo is 



e fe = lim E k = ^rA(k 

F^oc 24 



1 ^ 2 _(fe_3)(fc+l) ) 



(16) 



where A(n) is the number of alternating sign matrices 

A (n) = Uk=a( Sk + l ) [ /( n + k ) [ - 9 J- From Stirling's ex- 
pansion, the asymptotic behaviour of this expression is 
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